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The Littlest Higgs Model with T-parity is one of the attractive candidates of physics 
beyond the Standard Model. One of the important predictions of the model is the exis- 
tence of new heavy gauge bosons, where they acquire mass terms through the breaking 
of global symmetry necessarily imposed on the model. The determination of the masses 
are, hence, quite important to test the model. In this paper, the measurement accuracy 
of the heavy gauge bosons at ILC is reported. 

1 Introduction 

There are a number of scenarios for new physics beyond the Standard Model. The most 
famous one is the supersymmetric scenario. Recently, alternative one called the Little Higgs 
scenario has been proposed [TJ[2]. In this scenario, the Higgs boson is regarded as a pseudo 
Nambu-Goldstonc boson associated with a global symmetry at some higher scale. A 
symmetry called T-parity is imposed on the models to satisfy constraints from electroweak 
precision measurements [3J El [5] ■ Under the parity, new particles are assigned to be T-odd 
(i.e. with a T-parity of —1), while the SM particles arc T-even. The lightest T-odd particle 
is stable and provides a good candidate for dark matter. In this article, we focus on the 
Littlest Higgs model with T-parity as a simple and typical example of models implementing 
both the Little Higgs mechanism and T-parity. 

In order to test the Little Higgs model, precise determinations of properties of Little Higgs 
partners are mandatory, because these particles are directly related to the cancellation of 
quadratically divergent corrections to the Higgs mass term. In particular, measurements 
of heavy gauge boson masses, Little Higgs partners for gauge bosons, are quite important. 
Since heavy gauge bosons acquire mass terms through the breaking of the global symmetry, 
precise measurements of their masses allow us to determine the most important parameter 
of the model, namely the vacuum expectation value of the breaking. 

We studied the measurement accuracy of masses of the heavy gauge bosons at the inter- 
national linear collider (ILC). In addition, the sensitivity to the vacuum expectation value 
(f) was estimated. In this paper, the status of the study is shown, and the detail of this 
study is described in [6] . 



2 Representative point and target mode 

In order to perform a numerical simulation at ILC, we need to choose a representative point 
in the parameter space of the Littlest Higgs model with T-parity. Firstly, the model param- 
eters should satisfy the current electroweak precision data. In addition, the cosmological 
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e+e -> A H Z H 


e+e -> Z H Z H 


e+e" ■* W+PF H 


500 GcV 


1.91 (fb) 






1 TcV 


7.42 (fb) 


110 (fb) 


277 (fb) 



Table 1: Cross sections for the production of heavy gauge bosons. 



observation of dark matter relics also gives important information. Thus, we consider not 
only the electroweak precision measurements but also the WMAP observation [7] to choose 
a point in the parameter space. We have selected a representative point where Higgs mass 
and / arc 134 GeV and 580 GeV, respectively. At the representative point, we have obtained 
12dm'' 2 of 1.05. The masses of the heavy gauge bosons are (Ma h , M\y H , Mz H ) = (81.9 GeV, 
368 GeV, 369 GeV), where Ah, Zh, and Wh are the Little Higgs partners of a photon, Z 
boson, and W boson, respectively. Here, Ah plays the role of dark matter in this model 
[9] . Since all the heavy gauge bosons are lighter than 500 GeV, it is possible to generate 
them at ILC. 

There are four processes whose final states con- 
sist of two heavy gauge bosons: e+e" — > Ah Ah, ( a ) 
Ah^h, ZhZh, and W£. The first process is un- - 
detectable, thus not considered in this article. The 
cross sections of the other processes are shown in 
Table Q] Since mA H + mz H is less than 500 GeV, 
Ah^h can be produced at the ^fs = 500 GeV. At ^ 

t/s = 1 TeV, wc can observe WjJ"Wg with large cross g - w ~ g - _ 

section. We, hence, concentrate on e+e" — ► Ah-^h at \j' Z H ~T\» H 

yTs = 500 GcV and e+e" W+W^ at ^ = 1 ^ ^ + ^^JZ^^ 
TeV. Feynman diagrams for the signal processes are H 
shown in Fig. [TJ Note that Zn decays into A^_h, and 

Wi decays into AhW* with almost 100% branching Fi S ure ) : Diagrams for signal pro- 
fractions, cesses; (a) e+e -> A h ^h and (b) 




3 Simulation tools 



We have used MadGraph [10] to generate e+e" — * 

Ah^h at = 500 GeV, while e+e" — * W^W^ at -y/s = 1 TeV and all the standard model 
events have been generated by Physsim [TTj. We ignored the initial- and final-state radiation, 
beamstrahlung, and the beam energy spread for study of e+e - — > Ah^h at = 500 GeV, 
whereas their effects were considered for study of e+e~ — ► W^W-j^ at y's = 1 TeV where 
the beam energy spread is set to 0.14% for the electron beam and 0.07% for the positron 
beam. The finite crossing angle between the electron and positron beams was assumed to to 
be zero. In both event generators, the helicity amplitudes were calculated using the HELAS 
library [12] . which allows us to deal with the effect of gauge boson polarizations properly. 
Parton showering and hadronization have been carried out by using PYTHIA6.4 [13) . where 
final-state tau leptons are decayed by TAUOLA [14] in order to handle their polarizations 
correctly. The generated Monte Carlo events have been passed to a detector simulator called 
JSFQuickSimulator, which implements the GLD geometry and other detector-performance 
related parameters [15] . 
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Process 


Cross sec. [fb] 


# of events 


# of events after all cuts 


A H Z H -► A H A H bb 


1.05 


525 


272 


vvh — > i^fofo 


34.0 


17,000 


3,359 


Z/i — > wfofr 


5.57 


2,785 


1,406 




496 


248,000 


264 


ZZ — + 


25.5 


12,750 


178 


z/i/Z — > 


44.3 


22,150 


167 


7Z — > 766 


1,200 


600,000 


45 



Table 2: Signal and backgrounds processes considered in the AhZh analysis. 



4 Analysis 

In this section, we present simulation and analysis results for heavy gauge boson productions. 
The simulation has been performed at yfs = 500 GeV for the ArZh production and at s/s = 
1 TeV for the Wjj Wh production with an integrated luminosity of 500 fb -1 . 

4.1 e+e- -> AhZh at 500 GeV 

Ah and Zh are produced with the cross section of 1.9 fb 
at the center of mass energy of 500 GeV. Since Z-r de- 
cays into Ah and the Higgs boson, the signature is a 
single Higgs boson in the final state, mainly 2 jets from 
h — > bb (with a 55% branching ratio). We, therefore, de- 
fine AnZn — > AnAnbb as our signal event. For back- 
ground events, contribution from light quarks was not 
taken into account because such events can be rejected 
to negligible level after requiring the existence of two b- 
jets, assuming a ^-tagging efficiency of 80% for 6-jets with 
15% probability to misidentify a c-jet as a 6-jet. This b- 
tagging performance was estimated by the full simulation, 
assuming a typical ILC detector. Signal and background 
processes considered in this analysis are summarized in Figure 2: A typical event of 
Tabled Figure [2 shows a typical A H Z H event seen in the A H Z H in the simulator, 
detector simulator. 

The clusters in the calorimeters are combined to form a jet if the two clusters satisfy 
Uij < Vent- Uij is defined as 

_ 2E i E j (l-cos8 ij ) 

Vij " p2 ' V 1 -' 

vis 

where 6ij is the angle between momenta of two clusters, E^ are their energies, and E v { s is 
the total visible energy. All events are forced to have two jets by adjusting y cu t- We have 
selected events with the reconstructed Higgs mass in a window of 100 — 140 GeV. Since 
Higgs bosons coming from the WW fusion process have the transverse momentum (pt) 
mostly below W mass, pt is required to be above 80 GeV in order to suppress the vvh — > 
vvbb background. Finally, multiplying the efficiency of double 6-tagging (0.8x0.8 = 0.64), we 
are left with 272 signal and 5,419 background events as shown in Tabled which corresponds 
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E h [GeV] 



A H massfGeV) 



Figure 3: (a)Energy distribution of the reconstructed Higgs bosons with remaining back- 
grounds after the mass cut. (b) Probability contours corresponding to 1- and 2-cr deviations 
from the best fit point in the An and Zn mass-plane. The shaded area shows the unphysical 
region of m^ H + m ZH > 500 GeV. 

to a signal significance of 3.7 (= 272/V5419) standard deviations. The indication of the 
new physics signal can hence be obtained at yfs = 500 GeV. 

The masses of An and Zn bosons can be estimated from the edges of the distribution of 
the reconstructed Higgs boson energies. This is because the maximum and minimum Higgs 
boson energies {E max and E m i u ) arc written in terms of these masses, 

E m in = lz H E% - PzhIzhPh, (2) 

where (3z H (jz H ) is the ^(7) factor of the Zn boson in the laboratory frame, while is 
the energy (momentum) of the Higgs boson in the rest frame of the Z-r boson. Note that 
E* h is given as (Mf H + M* - M\ h )/(2M Zh ). 

Figure [3Ja) shows the energy distribution of the reconstructed Higgs bosons with re- 
maining backgrounds. The background events are subtracted from Fig. [3ja), assuming that 
the background distribution can be understand completely. Then, the endpoints, E ma , x and 
-E m in, have been estimated by fitting the distribution with a line shape determined by a 
high statistics signal sample. The fit resulted in itla h and mz H to be 83.2 ± 13.3 GeV and 
366.0 ± 16.0 GeV, respectively, which should be compared to their true values: 81.85 GeV 
and 368.2 GeV. Figure [3[b) shows the probability contours for the masses of Ah and Z-r. 

Since the masses of the heavy gauge bosons are from the vacuum expectation value (/) , 
/ can be determined by fitting the energy distribution of the reconstructed Higgs bosons. 
Then, / was determined to be / = 576.0 ± 25.0 GeV. 

4.2 e+e- -> W^W^ at 1 TeV 

W# Wft production has large cross section (277 fb) at ILC with ^fs = 1 TeV. Since W,f de- 
cays into Ah and with the 100% branching ratio, analysis procedure depends on the 
W decay modes. In this analysis, we have used 4- jet final states from hadronic decays of 
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Process 


cross sec. [fb] 


# of events 


# of events after all cuts 


W n W H ~* AnA K qqqq 


106.5 


53,258 


37,560 


W+W- -> gggg 


1773.5 


886,770 


306 


e + e~W + W~ — > e + e~qqqq 


464.9 


232,442 


23 


ev e WZ — > ev e qqqq 


25.5 


12,770 


3,696 


Z R Zn -> A H A H /i/i 


99.5 


49,757 


3,351 


vvW + W~ — > vvqqqq 


6.5 


3,227 


1,486 



Table 3: Signal and background processes considered in the Wjj W H analysis. 



two W bosons, WftWjH — > Ah 9999- Signal and background processes considered in the 
analysis are summarized in Table [3l 

All events have been reconstructed as 4- jet events by adjusting the cut on y-values. In 
order to identify the two W bosons from W^ decays, two jet-pairs have been selected so as 
to minimize a % function, 

X 2 = C C M W1 - ^M w f/a 2 Mw + C c Mw2 - tr M w ) 2 /a 2 Mw , (3) 

where rGC Mjyi(2) is the invariant mass of the first (second) 2-jet system paired as a W 
candidate, tr Mw is the true W mass (80.4 GeV), and (Jm w is the resolution for the W mass 
(4 GeV). We required \ 2 < 26 to obtain well-reconstructed events. Since An bosons escape 
from detection resulting in a missing momentum, the missing transverse momentum ( miss pr) 
of the signal peaks at around 175 GeV. We have thus selected events with mlss p T above 84 
GeV. Then, the reconstructed W energy is required to be between GeV to 500 GeV. The 
numbers of events after the selection cuts are shown in Table |3l The number of remaining 
background events is much smaller than that of the signal. 

As in the case of the AnZn production, the masses of An and Wg. bosons can be 
determined from the edges of the W energy distribution. Figure Ufa) shows the energy dis- 
tribution of the reconstructed W bosons. After subtracting the backgrounds from Fig[4ja), 
the distribution has been fitted with a line shape function. The fitted masses of An and 
Wh bosons are 82.29 ± 1.10 GeV and 367.8 ± 0.8 GeV, respectively, which are to be com- 
pared to their input values: 81.85 GeV and 368.2 GeV. Figure 0Jb) shows the probability 
contours for the masses of An and Wn at 1 TcV. The mass resolution improves dramatically 
at = 1 TcV, compared to that at y/s = 500 GeV. Then, / = 579.7 ± 1.1 GeV was 
obtained by fitting the energy distribution of the reconstructed W bosons. 

5 Summary 

The Littlest Higgs Model with T-parity is one of the attractive candidates of physics beyond 
the Standard Model since it solves both the little hierarchy and dark matter problems simul- 
taneously. One of the important predictions of the model is the existence of new heavy gauge 
bosons, where they acquire mass terms through the breaking of global symmetry necessarily 
imposed on the model. The determination of the masses are, hence, quite important to test 
the model. 

We have performed Monte Carlo simulations in order to estimate measurement accuracy 
of the masses of the heavy gauge bosons at ILC. At ILC with y/s = 500 GeV, it is possible 
to produce An and Zn bosons. Here, we can observe the excess by AnZn events in the 
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Figure 4: (a) The energy distribution of the reconstructed W bosons with remaining back- 
grounds after the selection cuts, (b) Probability contours corresponding to 1-, 3-, and 5-a 
deviations in the Ah and Wh mass-plane. 



Higgs energy distribution with the statistical significance of 3.7-sigma. Furthermore, the 
masses of these bosons can be determined with accuracies of 16.2% for Ah and 4.3% for Z^. 
Once ILC energy reaches y/s = 1 TeV, the process e + e~ — > Wjj W# opens. Since the cross 
section of the process is large, the masses of Wh and Ah can be determined as accurately 
as 1.3% and 0.2%, respectively. Then, the vacuum expectation value, /, can be determined 
with accuracy of 4.3% at y/s = 500 GcV and 0.2% at 1 TeV. 
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